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ABSTRACT 

Radio observations of the afterglow of the 7-ray burst GRB970508 provide unique 
new constraints on afterglow models. The quenching of diffractive scintillation at 
~ 4week delay provides the first direct estimate of source size and expansion rate. 
It implies an apparent size R ~ lO^^cm and expansion at a speed comparable 
to that of light at t ~ 4weeks, in agreement with the fireball model prediction 
R = 10^''(t/week)^/^cm. The radio flux and its dependence on time and frequency at 
1-5 week delay are in agreement with the model and imply a fireball energy (assuming 
spherical symmetry) ~ lO^'^erg, consistent with the value inferred from observations 
at shorter delay. The observed radio behavior deviates from model predictions at 
delays > 5weeks. This is expected, since at this delay the fireball is in transition from 
highly-relativistic to sub-relativistic expansion, with Lorentz factor 7 < 2. Deviation 
may be due to a change in the physical processes associated with the shock wave as it 
becomes sub-relativistic (e.g. a decrease in the fraction of energy carried by magnetic 
field), or to the fireball being a cone of opening angle ~ I/7 ~ 1/2. We predict the 
future behavior of the radio fiux assuming that the latter interpretation is valid. These 
predictions may be tested by radio observations in the frequency range 0.1~10GHz on 
time scale of months. 
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Introduction 



The availability of accurate positions for GRBs from the BeppoSAX satellite ( Posta et al. 



1997a| , Feroci et al. 1997 , Heise et al. 1997 , Costa et al. 1997c ) allowed for the first time to detect 



delayed emission associated with GRBs in X-ray ( posta et al. 1997a , b, Feroci et al. 1997 , Pirc 
let al. 1997ai b), optical (|Groot et al. 199% |Sahu et al. 1997 , [van Paradijs et al. 1997|, |Bond| 



19971 , pjorgovski et al. 1997| ) and radio ( [Frail fc Kulkarni 19971) wave-bands. The detection 



of absorption lines in the optical afterglow of GRB970508 provided the first direct estimate of 
source distance, constraining the redshift of GRB970508 to 0.86 < z < 2.3 ( [Metzger et al. 1997 ) 
Observed X-ray to radio afterglow are most naturally explained by models based on relativistic 
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blast-waves at cosmological distances (Paczynski & Rhoads 1993, Meszaros & Rees 1997, Vietri 



19974 [Waxman 19974 b, [Wijers, Rees &: Meszaros 1997| ). Using these models, combined radio and 
optical data allowed for the first time to directly estimate the total GRB energy, implying an 
energy of ~ 10^^ erg (assuming spherical symmetry) for GRB970508 ( Waxman 1997b| ). 

Radio observations of the afterglow of GRB970508 ( [Frail et al. 1997| , [Taylor et al. 1997| ) 
provide unique new information on the afterglow source. Shortly after the first detection of GRB 



radio afterglow, it was pointed out by Goodman ( Goodman 1997 ) that if the source angular size 
is as small as predicted by fireball models, ~ 1/xas, then the radio flux should be modulated by 
scintillation due to the local inter-stellar medium. The predicted modulation has been observed 
( Frail et al. 1997 , Taylor et al. 1997| ) and it provides for the first time direct constraints on 
source size and expansion rate. We discuss in §3 the implications to fireball models of the 
observed modulation. In §4 we discuss the implications of the information on long term afterglow 
behavior provided by the radio monitoring of GRB970508 over three months. Deviation from 
model predictions is expected on this time scale ( [Waxman 1997b| ), as the fireball decelerates from 
highly-relativistic to marginally relativistic speed. However, the lack of a theory describing the 
relevant physical processes associated with the shock wave (magnetic field generation, energy 
transfer to electrons) does not allow a unique interpretation of the observed behavior. We 
derive equations describing the fireball evolution in the non-relativistic regime, which allow 
to demonstrate that the observed behavior can not be accounted for by the deviation of the 
hydrodynamic evolution from that predicted by the highly relativistic scalings. We discuss possible 
interpretations of the long term behavior, and predict the future fireball radio emission under 
the assumption that the observed behavior is due to a finite opening angle of the fireball. We 
summarize our conclusions and predictions in §5. 

We note here that two types of fireball models were considered in the literature as possible 



interpretations of existing GRB970508 afterglow data. In the "adiabatic" model (Waxman 1997b) 
the fireball radiates over time only a small fraction of its energy, which is of order lO^^erg. In the 
"radiative" models ( Vietri 1997b[ , Katz & Piran 1997) the fireball radiates most of its energy over 
a short time, and its energy is reduced to ^ 10^^ erg on time scale of days. We discuss in this 
paper mainly implications of radio observations to the adiabatic model since, as we demonstrate 
in §4, long term radio observations rule out the radiative models, and since, as we briefiy show in 
52 below, the radiative models are not self-consistent. 



2. Radiative versus non-radiative models 

In fireball afterglow models, a highly relativistic shell encounters, after producing the GRB, 
some external medium. As the shell decelerates it drives a relativistic shock into the surrounding 
medium. This shock continuously heats fresh gas and accelerates relativistic electrons, which 
produce the observed radiation through synchrotron emission. The magnetic field B behind the 
shock and the characteristic Lorentz factor 7e of the electrons are determined in these models by 
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assuming that fractions and of the fireball energy are carried by electrons and magnetic field 
respectively. This implies S^/Svr = A^Bl'^nmpc'^ and 7e = 6.e7^p/n^e (here n is the number density 
of the surrounding medium, and 7 is the fireball Lorentz factor). There is no theory which allows 
to determine the values of and ^s- Afterglow observations are consistent with ~ '^B ~ 0.1 
during the relativistic expansion of the fireball ( |Waxman 1997b( ) . 

The underlying assumption of radiative models is that all the kinetic energy lost by the 
fireball as it decelerates is radiated away. This requires an efficient process converting the kinetic 
energy flux to electron thermal energy, and requires the electron cooling time to be much shorter 
than the dynamical time, i.e. the time scale for expansion in the fireball rest frame Td = r/A^yc 
(here r is the fireball radius). The ratio of synchrotron cooling time, = Girniec/aTjeB'^ , to 
dynamical time is 



2 



e = - = - — iUBCTTurY)-' = 16iUB/0.25r'inim7T, (1) 
Td 4 yrup J 

where n = Inicm"^, r = 10^''ri7cm. (This result is similar to that given in jVietri 1997bi ). In 
the radiative model of Katz & Piran (Katz & Piran 1997) the fireball expands into a uniform 
density medium, ni = 1, and due to rapid energy loss the fireball becomes non relativistic at 6 
day delay with riy ~ 7 ~ 1. In the radiative model of Vietri ( Vietri 1997b), the ambient medium 
density drops with radius as r~'^ allowing the fireball to remain relativistic as it loses energy. In 
this case ni 10^^ and ru ~ 7 3 at t = 6day. In both models, therefore, 1 at this time. 
Furthermore, the time dependence of G, oc for n oc r~'^ and Q oc t^^^ for uniform density, 
implies that the assumption G ^ 1 is not valid in both models for t ^ Ihr. 



3. Source size and expansion rate 

Due to relativistic beaming, the radiation from a relativistic fireball seen by a distant observer 
is emitted from a cone of the fireball around the source-observer line of sight, with an opening 
angle ~ I/7. The apparent radius of the emitting cone is R = r/7, where r is the fireball radius, 
and photons emitted from such a cone are delayed, compared to those emitted on the line of sight, 
by t = r/27^c. Thus, the apparent radius of the fireball is R = 2j{r)ct (a detailed calculation of 



fireball emission [ Waxman 1997c introduces only a small correction, R = 1.9j{r)ct), where r and 



t are related by t = r/27^c. Using eqs. (1) & (2) of ( [Waxman 1997b| ) we have 
and 

•1 + A'/' fE^Y 



i? = 8 X lO^^ j (^j (3) 

Here E = 10^^i^52erg is the fireball energy and t = It^week. (Note, that the relation t = r/2j'^c 
should not be replaced hy t = r/167^c, as recently argued in [pari 1997 |, since the latter relation 
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holds only for the arrival time of photons which are emitted on the line of sight and not valid for 
most of the photons, which are emitted from a cone of opening ang le ~ 1/7 [IWaxman 1997^ ). 



Scattering by irregularities in the local interstellar medium (ISM) may modulate the 
observed fireball radio flux ( Poodman 1997 ). If scattering produces multiple images of the 



source, interference between the multiple images may produce a diffraction pattern (on an 
imaginary plane perpendicular to the line of sight), leading to strong variations of the flux as 
the observer moves through the pattern. Two conditions need to be met in order to produce 
such diffractive scintillation: (i) The scattering should be strong enough to produce multiple 
images; (ii) The source size should be small enough so that different points on the source produce 
similar diffraction patterns. There is significant observational evidence that ISM electron density 
fluctuations are described by a power-law spectrum, {6Ne{k)6Ne{—k)) = C'^k~^^/^, where k is 
the spatial wave number. For this distribution, the characteristic deflection angle is given by 
9d = 2.34A^^/^re^^(<S'M)^/^, where A is the wavelength, the classical electron radius, and the 
scattering measure SM is the integral of along the line of sight. In the frame work of the "thin 
screen" approximations, i.e. assuming that all scattering occurs in a narrow layer at distance dgc, 
multiple images are produced for frequencies ( Goodman 1997| ) 



''<Mi'H '!L,. )'"'ghz. (4) 



sckpc Vl0-3.5m-20/3kpc 

Here, values were chosen for the ISM scattering properties SM and dsc, which are typical for 
sources at high Galactic latitude (6 = 27° for GRB970508). The characteristic length scale of the 
diffraction pattern is 

A _ .^10 6/5/ SM 



dx = — — = 3.3 X lO'^'u"'" — cm. (5) 

In order for the diffraction patterns produced by different points on the source to be similar, so 
that the pattern is not smoothed out due to large source size, the angular source size Og should 
satisfy Ogdsc < dx. For a source at redshift z = 1 this requirement implies an upper limit to the 
apparent source size 

where = lOi^ioGHz, hj^ is the Hubble constant in units of 75km/sMpc. Due to the weak 
dependence of the angular diameter distance on z, the upper limit on R is not sensitive to source 
redshift. The main uncertainty in (|6|) is due to uncertainty in the scattering properties of the ISM. 
Combining (P) and ^ we flnd that 

/ S'/Vf \ ^3/17 

i? < 1.1 X lO^^d^^Z'h^i ( io„3.5^-20/3,pJ 

is required to allow diffractive scintillation. Due to variations in ISM scattering properties along 
different lines of sight, the numerical value in (^) is accurate to a factor of a few. 
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The frequency range Az/ over which the diffraction pattern is similar, and therefore over 
which flux modulation is correlated, is 

C _ ( SM \-6/5 / y \22/5 

= = °-<'-P- l l0-3.5n,-20/3i,p J GHz. (8) 

For a characteristic velocity through the diffraction pattern of ~ 30km s^^, due to Earth's orbital 
motion and to the Sun's peculiar velocity, implies a time scale for variations ~ 3hr. 

Comparing (|3|) and (^) we find that on time scale of weeks the apparent fireball size is 
comparable to the maximum size for which diffractive scintillation is possible. On shorter time 
scale, therefore, strong modulation of the radio flux is expected. On longer time scale we expect 
diffractive scintillation to be quenched due to large source size. When diffractive scintillation 
is quenched, the flux is nevertheless expected to be modulated due to refraction (i.e. due to 
focusing/defocusing of rays). However, the modulation amplitude should decrease (to ~ 10%, 



Goodman 1997 ), and modulation should be correlated over a wide frequency range. Figure 1 
presents the light curves of the radio afterglow at 8.46GHz, 4.86GHz and 1.43GHz. Figure 2 
presents the fluxes at 4.86GHz as a function of the flux at 8.46GHz. Strong modulation of the 
radio flux is observed during the first month, accompanied by strong variations in the ratio of 
flux at 8.46GHz and 4.86GHz. This behavior is consistent with that expected due to diffractive 
scintillation. The radio flux is not sampled at a sufficient rate to determine whether the variability 
time scale is consistent with that expected for diffractive scintillation. At delays longer than ~ 1 
month the modulation amplitude decreases, and the flux ratio is consistent with being constant. 
This is consistent with quenching of diffractive scintillation due to increased source size. 

Observations are therefore in agreement with fireball model predictions. They imply that the 
source size is close to the upper limit given by (0) after ~ Imonth. This is consistent with (^) 
and implies expansion at a speed comparable to that of light. Due to the very weak dependence 
of R on fireball model parameters it is not possible to accurately determine parameters based on 
the quenching of diffractive scintillation. On the other-hand, since R is very insensitive to model 
parameters, reducing the uncertainty in the size estimate based on scintillation, by reducing the 
uncertainty in ISM scattering properties towards the GRB, would provide a stringent test of the 
fireball model. 



4. Long term behavior 
4.1. Relativistic regime 



The radio light curves are compared with the predictions of the adiabatic fireball model 
( Waxman 1997b| ) in Figure 1. In this model, the observed frequency at which the synchrotron 
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spectral intensity peaks is 



.R 



5 X 10 



12 



1 + Z 



1/2 



and the observed intensity at is 



R 



l + z 



1-1/V2 



n 



1/2 



(CB/0.1)i/2i?52mJy. 



(9) 



(10) 



The superscript R imphes that the expressions are vahd for the highly relativistic regime. The 
flux at is produced by electrons at the characteristic electron energy, e = jeirieC^- The flux 
at higher frequency is produced by higher energy electrons. For a power-law electron spectrum, 
dNe/dSe oc e~P, Fi, oc at > Vm- Typical parameters required to fit observations are 

E52 ~ f^i ~ 1, 'Ce ~ ~ 0.1, and p ~ 2. 

The flux at low frequency, < fm, is due to the extension of synchrotron emission of electrons 
at e = 7emeC^ to frequencies v < Vm [Fy oc {v /vmY^"^ a± v ^ Vm]- At low frequency self-absorption 
becomes significant. The self-absorption frequency, where the fireball optical depth is unity, is 



1 



l + z 



(ee/0.2)-i(eB/0.1)i/5i?5f nf^GHz, 



(11) 



and at V < Urn 
a factor (1 



the optical depth is given by = {yjvA) 



-5/3 



Self-absorption reduces the flux by 



The solid smooth curves in Figure 1 give the model fluxes for £'52 = ni = 2, = 0.2 and 



= 0.1. These are essentially the same parameter values inferred in (Waxman 1997b) from 
optical and radio afterglow data at delays t < 6day. Although many simplifying assumptions were 
made, in order to obtain a simple description of fireball behavior, the model is in agreement with 
the data obtained during the first ~ 5 weeks. At later time model curves deviate from the data. 
This behavior is expected ( [Waxman 1997b ), as the fireball decelerates from highly-relativistic to 
marginally relativistic speed. At this stage the scalings (§J3|) do not give an accurate description of 
the fireball dynamics. Furthermore, there is no theory which allows to determine the parameters 

and S,B- These parameters may change as the shock decelerates, thus affecting the predictions 
(^-|ll). Prior to discussing in §5 the possible implications of observations at t > Sweeks we derive 
in §4.2 the equations describing the fireball dynamics at the non-relativistic stage. This would 
allow to estimate the effects due to deviation from the scaling laws (P,^). 



4.2. Transition to the non-relativistic regime 



As the fireball becomes non-relativistic its expansion approaches that described by the 
Sedov-von Neumann- Taylor solutions ( |Sedov 1946| , [von Neumann 1947 , p?aylor 1950 ). At this stage 
the shock radius is given by r = S^Q{'y){Et'^ /nnipC^y^^ , where ^0 is a function of the adiabatic index 
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of the gas 7. .^0 = 0.99 for 7 = 4/3 (relativistic fluid) and ^0 = 1-15 for 7 = 5/3 (non-relativistic 
limit). The non-relativistic behavior may be described as 

/3 = r/c = S\T/rMR)-^'\ r/vMR = (eo/1.15)(^/^iVi^)'/^ (12) 

by defining 

r^R = 1.0 X W'^ (^^^ cm, t^R = 2i^-^ (^^^ week. (13) 

With these definitions, the time dependence (||) of the Lorentz factor in the highly-relativistic 
regime may be written as 7 = 1.3{t/tNR)~^^^ ■ Thus, the relativistic solution is valid for t <^ t^R, 
where 7^^ ^ 1, and the non-relativistic solution for t ^ t^ji, where <^ 1 



For the non-relativistic regime, assuming that fractions and of the dissipated energy 
are carried by magnetic field and electrons imply B'^/Sn = pi^nrUpC? and characteristic electron 
Lorentz factor 7e = (mp/2me)^e/3^- The frequency at which the synchrotron intensity peaks, 
= 7eei3/27rmeC, is 

,NR _ 



^ + \ee/0.2)2(e5/0.1)V2ni/35GHz. (14) 



The intensity at is F^J^ = Ne{V^e^B /2'KmeC^){l + z)/A7rdl, where Ne = Anr^n/S is the 
number of radiating electrons, dL the luminosity distance. Using (QOh we have 



.rNR. 

Extrapolation of the non-relativistic expression (15) to t = tNR gives a peak flux similar to that 



given by (10) for the relativistic regime. This is expected, since in the relativistic regime F^^ 
is independent of time, and therefore of 7. Thus, as the fireball decelerates to non-relativistic 
speed the peak flux is approximately given by the relativistic expression (0). At lat er time, 
t ^ tNR, Fy^ oc (t/tNR)^^^ ■ Extrapolation of the relativistic expression (P) for i^m to t = tjsfR 
gives Um = 100(^e/0.2)^(^5/0.1)^/^n}''^GHz, significantly higher than the extrapolation of the 
non-relativistic expression (p!4|). Thus, as the fireball decelerates Vm decreases with time faster 
than given by (|9|), and Um oc (t/t7Vi?)~^ for t S> t]\fR- 

The self-absorption frequency, where the fireball optical depth is unity, is 

.r = ^5r'/'(^Y'. (16) 



Comparing (|16D and (11) we find that the self-absorption frequency is time independent and is 
given by the relativistic expression for t < tATR. At later time ([T6| ) implies that the self-absorption 
frequency increases. However, ( [T6| ) is valid only as long as ua < ^m- Since the dependence of Vm 
on time for t > t^R is stronger than that of va, drops below va when va is not significantly 
higher than the value given by (0). At later time va decreases with time. For an electron 
spectrum dNe/dSe oc e~^, va oc {t/iNR)'"^^^ ■ 
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5. Implications 
5.1. Ruling out radiative expansion 



From the analysis of §4.2, the relativistic expression ( |10D is a good approximation for the 

fireball peak flux not only for 7 » 1, but as long as /3 ~ 1. Since radio observations imply 

that the fireball expands with /3 ~ 1 on time scale of weeks, the observed flux of order ImJy 
1/2 

implies £52^1 ^ 1 on weeks times scale. This rules out the radiative models, where the fireball 



energy decreases to lO'^^erg on day time scale. We note here that it was argued in ( [Katz fc Piran 



1997) that a mildly relativistic fireball, 7 — 1 ~ 1, with E ~ 10^^ erg would produce the observed 



O.lmJy flux at 1.43GHz at ~ 6day delay. This is in contradiction with our results, (10) and 



(JTsD, which imply a much lower flux. The ( |Katz &: Piran 1997 ) derivation is, however, not self 



consistent. The flux at 1.43GHz is obtained in ( Katz &: Piran 19971) using the Rayleigh-Jeans 
law, for which Fi, oc v'^ and which is valid only for high optical depth. From ([Tl|) and (|l^) the 



self-absorption frequency for the parameters chosen in ( Katz Sz Piran 1997| ), E ~ 10^^ erg and 



ni = 1, is ~ 0.2GHz. Using the Rayleigh-Jeans formula for the flux at 1.43GHz, where ~ 0.05, 
overestimates the flux by a factor (1 — e~'^")^^ ~ 30. 



5.2. Deviations during the transition to non relativistic expansion 

On time scale ^ Sweeks the observed radio behavior deviates from model predictions. The 
flux at 4.86GHz and 8.46GHz is significantly below model predictions at t ~ lOweeks. The 
increase of flux at 1.43GHz at this time indicates that the frequency at which the intensity 
peaks drops to ~ 5GHz at t ~ lOweek, where (^) predicts Um ~ lOOGHz. From the analysis 
of the previous section, this behavior can not be explained based only on the deviation from 
the highly-relativistic scaling laws (0J3|) as the fireball decelerates to mildly relativistic velocity, 



7 — 1 ~ 1. The synchrotron peak intensity is not expected to decrease [cf. eq. (15)], and therefore 
the decrease in 4.86GHz and 8.46GHz flux can not be accounted for. The peak frequency fm is 
expected to decrease faster than predicted by (P). However, it is expected to decrease to ~ 5GHz 
only when the fireball becomes non-relativistic, i.e. at t ~ t^ji ~ 24week [cf. eqs. ([l^, (|l2|)]. 
Note, that a change in the ambient medium density n is not likely to account for the observed 
behavior. A decrease in flux may result from a decrease in n (|lO|). However, the peak frequency 
@ is independent of n, and the apparent decrease in can not be accounted for. 

Clearly, the observed behavior may be explained by deviations of the equipartition fractions 
and from the values ~ ~ 0.1, which are implied by observations at t < 5weeks. Due to 
the lack of a theory determining these parameters, it is not possible to predict their dependence on 
shock Lorentz factor. Thus, if the observed behavior is due to changes in and ^s, it is difficult 
to predict the future fireball behavior. However, the observed deviations from the model may also 
result from a different effect. This is discussed below. 
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5.3. Non-spherical fireballs 

We have so far assumed that the fireball is spherically symmetric. The results are valid also 
for the case where the fireball is a cone of finite opening angle 9, as long as 7 > 1/^. In this case, 
the fireball energy E in the equations should be understood as the energy the fireball would have 
had if it were spherically symmetric. The actual fireball energy is E' ~ 9'^E/2. Deviations from 
the spherical model would appear at late time, as 7 decreases below 1/6 and the fireball starts 
expanding transversely as well as radially. Let us briefiy discuss the expected behavior at later 
time. 

After a transition phase the fireball would approach spherically symmetric expansion, which 
may again be described by the equations derived above, with E replaced by the actual fireball 
energy E' . Numerical calculations would probably be required to describe the fireball evolution 
in the stage it approaches spherically symmetric behavior. Qualitatively, as the fireball expands 
transversely its energy per unit solid angle along the line of sight decreases. Thus, it would appear 
as if the fireball energy is decreasing with time. This would lead to decrease in the peak flux 
[cf. ([lO|) ] and in the peak frequency [cf. @], in qualitative agreement with the observed trends. 
The time scale for transition to spherically symmetric behavior may be estimated as follows. The 
scaling laws (||) are derived from energy conservation, E ~ [Airr^ /3)nmpc'^ . This implies that 
at the radius rg, where 7 = 1/6, the rest mass energy contained in the sphere with r = rg is 
comparable to the total fireball energy E' . Thus, as the fireball approaches spherically symmetric 
behavior it necessarily becomes sub-relativistic, and at later time it is described by (1^16), with 



E replaced by E' . The transition to spherical non-relativistic behavior occurs on a time scale [cf. 
eq. (0)] tNR = 18(e2^52/ni)i/3^eek. 

The deviation at t ~ 5week from model predictions may therefore be accounted for by the 
fireball being a cone of opening angle 6 ~ l/7(5week) ~ 1/2. This implies that the "real" fireball 
energy is £" ~ 2 X lO^^erg. In this case, the decrease in flux and in peak frequency are accounted 
for by the transition to spherical non-relativistic behavior. The transition should take place on 
a time scale t^R ~ 12week. On this time scale the peak frequency is expected to decrease to 
~ 4GHz [cf. eq. and the peak intensity to ~ O.SmJy [cf. eq. (|l5|)]. This is in agreement 

with the observed behavior. The fireball behavior for t S> 12week is described by (|l^-p!^), with 
E52 ^ 0.2, ni = 2. 



6. Conclusions 

Comparison of radio observations of the afterglow of the 7-ray burst GRB970508 with fireball 
model predictions lead to the following conclusions: 

• The source size implied by the quenching of diffractive scintillation at ~ 4week, ~ lO^^cm, 
and the inferred expansion at a speed comparable to that of light are consistent with 
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published ( [Waxman 1997b| ) model predictions @, (|3|). 

The radio flux and its dependence on time and frequency at 1-5 week delay are in agreement 
with the model (Figure 1) and imply a fireball energy (assuming spherical symmetry) 
~ lO^^erg. This is consistent with the value inferred from observations at shorter delay 
( Waxman 1997b ), and rules out "radiative" models ([Vietri 1997b , Katz &: Piran 1997 ), 



where fireball energy is reduced to ~ 10 erg on day time scale. 

• The deviation of observed radio behavior from model predictions at delays > 4weeks 
is expected, as on this time scale the fireball is in transition from highly-relativistic to 
sub-relativistic expansion, with Lorentz factor 7 < 2. We have shown that the observed 
behavior can not be accounted for by the deviation of the hydrodynamic behavior from that 
predicted by the highly relativistic scalings (|2|,^), or by changes in ambient medium density. 

• The observed behavior may be explained by deviation of the equipartition fractions and 
S^B from the values ~ ~ 0.1, which are implied by observations at t < Sweeks. Due 
to the lack of a theory determining these parameters, it is not possible to predict their 
dependence on shock Lorentz factor. Thus, if the observed behavior is due to changes in 
and ^B, it is difficult to predict the future fireball behavior. 

• However, the observed behavior at t > 5week may also be accounted for by the fireball being 
a cone of finite opening angle, 6 ~ l/7(5week) ~ 1/2, which implies that the "real" fireball 
energy is E' ~ 6'^E/2 ~ 2 x lO^^erg. In this case, at t ~ 5week the fireball rapidly expands 
transversely, leading to a decrease in the energy per solid angle along the line of sight. This 
may account for the observed decrease in peak flux and peak frequency. On time scale 
tNR ~ 12week [cf. eq. (|l3|)] the flreball approaches spherical non-relativistic expansion. On 
this time scale the peak frequency is expected to decrease to ~ 4GHz [cf. eq. (p^)], and the 
peak intensity to ~ O.SmJy [cf. eq. (^)]. This is in agreement with the observed behavior. 

If the observed behavior is indeed due to the flreball being a cone of flnite opening angle, than 
the future behavior may be predicted. On time scale of ~ 12week, the flreball should approach 
spherical non-relativistic expansion. At times t ^ 12week the behavior is given by eqs. (|l2|-|l^), 



which describe the spherical non-relativistic behavior, with £^52 — 0.2, ni = 2. The frequency at 
which the intensity peaks should decrease with time, oc i^^, dropping to 0.3GHz in ~ 0.5yr. 
The peak intensity should be ~ O.SmJy over this period. If the electron energy distribution is 
similar to the distribution inferred for the relativistic regime, dNg/de oc with p ~ 2, then the 
flux aX. V > Vrn should decrease approximately as t~^, similar to the highly relativistic case, and at 
a given time the intensity should drop with frequency as z^~" with a ~ 0.5. 

The time dependence of the flux at 1/ > f m is similar in the non-relativistic and the relativistic 
regime. At early time, t < Iweek, the optical flux decrease as ( pjorgovski et al. 1997 ). 



However, the optical flux should drop at ~ 5week below the t ^ extrapolation from early time, as 
the effective flreball energy decreases. The flux expected at ~ lOweek in this model is m/j ~ 25. 



Agreement of the observed optical behavior with the behavior described above would provide 
support to the model discussed here, where it is assumed that both optical and radio fluxes are 
produced by the expanding fireball. 
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Fig. la. — Light curves of the radio afterglow of GRB970508 at 4.86GHz and 1.43GHz, compared 
with the predictions of the adiabatic fireball model ( |Waxman 1997b| ). 
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Fig. lb. — Light curve of the radio afterglow of GRB970508 at 8.46GHz, compared with firebah 
model predictions ( [Waxman 1997b ). 
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Fig. 2.— The flux at 4.86GHz as a function of the flux at 8.46GHz, for t > 43day and t < 43day 

(flrst/second half of observing time). The average value of the flux ratio / = 5'8.46/5'4.86 for 
t > 43day is / = 1.17 and the standard deviation in / is cr = 0.15 (the region of 1 standard 
deviation is bordered by the solid lines). The scatter is consistent with that expected due to errors 
in flux measurement (x^ = 17.6 for the hypothesis / = / for the 13 data points), and / is consistent 
with the value expected in the fireball model, / = (8.46/4.86)^/^ = 1.20. At delays t < 43day the 
average flux ratio is similar, / = 1.18. The scatter is, however, large, a = 0.45, and inconsistent 
with that expected from measurement errors (x^ = 267 for the hypothesis / = / for the 18 data 
points) . 



